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We explore the diffusion process in the non-Markovian spatio-temporal noise. There is a non- 
trivial short memory regime, i.e., the Markovian limit characterized by a scaling relation between 
the spatial and temporal correlation lengths. In this regime, a Fokker-Planck equation is derived by 
expanding the trajectory around the systematic motion and the non-Markovian nature amounts to 
. the systematic reduction of the potential. For a system with the potential barrier, this fact leads to 

' the renormalization of both the barrier height and collisional prefactor in the Kramers escape rate, 

Lj^ \ with the resultant rate showing a maximum at some scaling limit. 
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It has been a long-standing problem to systematically treat the spatio-temporal colored noise in the thermal 
diffusion^], 0, ||. This is in contrast to the successful development of the noise- induced diffusion processes in the 
cases of additive noisejl], multiplicative noiseQ, and temporally colored noise [f|. In the context of path coalescence, 
DeutschQ first examined the motion of a damped particle subjected to a force fluctuating in both space and time. 
Wilkinson, Mehlig and coworkers @, H, pursued a similar subject, intensively explored the generalized Ornstein- 
Uhlenbeck process in momentum space, and showed the anomalous diffusion where the spatial memory led to the 
staggered ladder spectraQ ■ A typical example of the spatio-temporal noise is given by the turbulent flow or a randomly 
moving gas which suspends small particles, and space coordinates often stand for collective or reaction coordinates 
and order parameters in general. Then small particles or molecules suspended by the turbulent flow can experience 
systematic forces. The above works 0,0, Hill) however, are limited to the diffusion with no systematic force. It is thus 
a natural challenge to consider the diffusion with a spatio-temporal correlated noise in the presence of the systematic 
force induced by a potential. The motivation of the present paper is to elucidate a scenario how the spatio-temporal 
correlated noise may radically affect the potential-induced motion with the purely-temporal noise. In this case, we 
shall see that there is a non-trivial Markovian limit characterized by the spatial and temporal correlation lengths, and 
the Fokker-Planck equation for the distribution of position, collective coordinates or the order parameter in general, 
is derived from the conservation of population. Then we shall investigate the Kramers escape rate problem where the 
potential renormalization plays a crucial role. The new escape rate matches up to renewed attentions paid for the 
Kramers escape rate theory [Tol ITT1. [l2j . 
ly-j \ This Letter is organized as follows: After a sketch of our model, the Markovian limit is discussed and then the 
. Fokker-Planck equation is derived from the conservation of the population. The escape rate is derived as an inverse 

of the mean first passage time. 
, Our model is described by the overdamped Langevin equation with a colored noise 

^ : V x(t) = -U'(x(t)) + f(x(£),t) 

a>: (f(x,t)) = o,(f(x,t)f(x',t')) = c(x-x',t-t'), (i) 

where 77 stands for the friction coefficient and the potential U (x) consists of a metastable well and a potential barrier 
(see Fig.l). The bracket () stands for the average over the noise f(x, t) and the initial distribution. The noise f(x, t) is 
assumed as a Gaussian process in accordance with the central limit theorem both in the spatial and temporal variables. 
An explicit expression of the noise f{x,t) is given in the later paragraph concerning our stochastic simulation. As a 
typical form of the autocorrelation function, we consider the Gaussian memory 

C(x,t) = Coe~$-& , (2) 

where £ is the spatial correlation length and r is the correlation time, respectively. The constant Co contains r and 
£ in general. Its asymptotic form is fully determined in a short correlation limit which we shall discuss below. We 



'Electronic address: monnai@suou.waseda.jp 
^Electronic address: sugita@a-phys.eng.osaka-cu.ac.jp 
t Electronic address: nakamura@a-phy.eng.osaka-cu.ac.jp 



2 



U(x) 



\ 

\ ^min 


^max 







X 



FIG. 1: The schematic illustration of the potential U(x). The thick line indicates the domain of the metastable state. 



note that essentially the same result is obtained for the other form of the memory, C( x ,t) = C s -^^e- t2 / T \ which 
behaves as Gaussian at the origin, and converges to the Dirac-delta. We assume that the usual overdamped Langevin 
equation is recovered as the position-independent noise case by taking the limit of £ — > oo and the temporal Markovian 
limit t — ► with Co = 2riksT /t^/tt. There is still a possibility of £ dependence as a Laurent series 

with the strongest divergence C^/^ n , (n > 0) where /£" is dimensionless and = 1. It will become clear 
that all the terms should vanish C^ n ' — for n > 0, by requiring that the diffusion coefficient is finite in the short 
correlation limit £ — > 0. At the moment, however, we employ the general expression ©.Note that Sinai's random walk 
in random environment leads to the anomalously slow diffusion [l3|, while in our case symmetric nature of random 
force may enhance the diffusion over the barrier. 

The Fokker-Planck equation at a given instant is obtained from the conservation of probability The stochastic 
Liouville equation for the probability distribution p{x 1 t) of each realization of noise is given as 

d Id 

T-p(x,t) = -Lp(x,t) —f(x,t)p(x,t), (4) 

at r\ ox 

where Lp(x,t) = —-j^-U' (x)p(x,t) is the deterministic Liouville operator. With use of the identity 



p(x,t) = e- tL p(x,0) - - I dse^- t ^^-f(x,s)p(x,s), (5) 



V Jo dx 



the equation (j4]) is rewritten in the form of diffusion equation as 



& d f 

= -Lp(x,t) - —-f(x,t)e~ tL p(x,0) 
at ox rj 



x,s)p(x,s). (6) 



By averaging over the noise, the population distribution P(x, t) = (p(x,t)) obeys the diffusion equation, 
Wt P ^ = Wx U '^ P ^ + ^ll 



ds(f(x,t)^f(y(x),s)p(y(x),s)). (7) 

Here, y(x) = e^^ s ^ L x — x + f s ds'U'(x(s'))/T] is the time-reversed deterministic trajectory which starts at x. The 
last term in Eq. (JT]) is evaluated as follows. After the averaging, we shall take the Markovian limit r — > 0, and 
thus the correlation between explicit noise terms and those at earlier times implicitly contained in the distribution 
p(y(x),s) is neglected. More precisely, the correlation between / and p is calculated by substituting the identity ([5]) 
into the stochastic Liouville equation twice, and applying the pair-wise factorization for Gaussian variables: 

'ds [ S ds> f ds"l/T 2 e-( t - s 'f/T 2 e-(°-s"f/r 2 
Jo Jo 

0(r)(r - 0), (8) 
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which turns out to be negligible in the Markovian limit. On the other hand, the correlation between noise terms is 

d S (f(x,t)(^f(y(x),s))) 
1 dy d _ (*-«< f )) 2 -I^jli. 



dsc (i + o(r)) 2 J^l(t - s)e -^+^-°r 



This means a systematic renormalization of the potential. We approximated the trajectory as y(x) = x + 
f s ds'U'(x(s'))/r) w X + U'(x)(t — s)/?7, because it appears in the exponentially decaying factor as t — s gets 
larger and, in the Markovian limit, only the short time integral does contribute to the correlation function. Choosing 
the constant Co as Eq. ((3|) so that the usual Langevin equation is recovered as a special case, the dominant contribution 
of Eq. ([9]) to the Fokker-Planck equation is 

u>(x) 2C^ V k B T 

V^(t^U'(x) 2 + ' 

where a Markovian limit is well-defined: Keeping the dimcnsionless parameter 

n n = 2C^rk B T/V^C +2 V (11) 

as a constant, both the temporal and spatial correlations go to zero, r, £ — * 0. Similarly the diffusion coefficient is 
given as 



k B T CW £77 



(12) 



which diverges in the limit t, £ — > for n > 0. Thus hereafter we consider the case of n = 0, where the diffusion 
coefficient approaches in the Markovian limit specified by k = «o- It is instructive to derive the diffusion constant 
in a more general way. We assume that the kernel C(x,t) rapidly vanishes as t/r —* 00. Then the time-integral of 
noise correlation function is calculated as 



t 

dsC{x — y(x), t — s) 



1 

V 2 Jo 

— / dsC{£ ,t ) 

v Jo £ T 

1 [* „„ U'{x)T t-s t-s 

— / dsC{S,— ,r ) 

V Jo 0/ r 

1 t-s. 



= ^, (13) 
ry 

where we employed the same approximation, i.e., the expansion around the deterministic trajectory y(x), and used 
the relation j = £ • y/irr]K/2kBT in the Markovian limit. The integral gives the diffusion coefficient. 
In the above-mentioned short correlation regime, the Fokker-Planck equation becomes 

|p(x, t) = i|- (u'(x){l -k) + k B T^) P(x, t). (14) 

It is one of our discoveries that the spatial correlation amounts to the systematic reduction of the drift velocity. The 
usual Langevin dynamics driven by the white noise is recovered in the case of long-enough spatial correlation length 
which is expressed as k — 0. Mathematically, this condition is satisfied when the ratio t/£ 2 — > 0, which only requires 
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that r should vanish faster than £ 2 . Intuitively, however, the spatial correlation length can be seen as infinite when it 
exceeds the typical length of the particles' displacement within the relaxation time r. The steady state distribution 
is now given as 

P S \x) = Ne -U(x)(l- K )/k B T^ p_gj 

which has the form of canonical distribution with the renormalized potential U(x)(l — k) (Fig. 3). Here, the natural 
boundary condition is assumed, and N is the normalization constant. 

Let us now proceed to the escape rate problem characterized by the potential in Fig.l. For the diffusion process 
governed by the Fokker-Planck equation (|14p . the mean first passage time is now straightforwardly calculated, and 
its inverse gives the escape rate. We assume the initial probability distribution P(x, t = 0) = S(x — xq) confined in a 
metastable region bordered by the free energy barrier at x = x max and the infinite^wall at x — — oo. The mean first 
passage time t(xo) from the initial position xq is given by the standard procedure [H as a quadratic integral 

r (a; ) = -H- [ Xma * dye u(y)(i- K )/k B T f V dze _cr W (i-«)/fc B r i (16) 

>ZbT J Xg J-oo 

where the dependence on the initial point may not be important at low temperature. For a high-enough barrier, the 
escape rate is evaluated by the steepest descent approximation at the maximum x ma x and the minimum x m in' 

r = r^o)- 1 = -^-^U"(x min )\U"(x max )\(l - K )e-^ (1 - K) . (17) 
2m] 

To confirm the validity of Eq. (|17jl , we numerically solved the non-Markovian Langevin equation {1} in the special 
case of Landau potential U = x 4 /4— x 2 /2, under sufficiently short correlations of noise. The spatio-temporal correlated 
noise can be given in a product form f(x, t) = gi{x,£)g2(t,r) where the noises 31=1,2(5, a) have the Gaussian correlation 

(gi(s,a)gj(s' , a)) oc b^e <? with the correlation a. In order to construct the colored noise gi[s,a), we considered 
an assembly of harmonic oscillators which mimics the Langevin force 

/ \ sinwfcS 
gi{s,a) = > (xkcoswks +pk ), (18) 

k 

where the coefficients Xk and pk are the Gaussian stochastic variables with zero means and the variances 

(x 2 k ) = di/mul 

(pl) = Oi/m, (19) 
with the mass m = 1 and the noise strength 8. The density of the states Gi(u>) is chosen as 

Gi(uj) tx Lj 2 e-^. (20) 

Then the desired Gaussian correlation with correlation length a is achieved. Note that the stochastic process f(x(i), t) 
is not uniquely determined from the variance, and the product form is one of the possible choices which guarantee the 
escape rate formula p7|) . We fixed the parameters r = 0.01, 77 = 1, fc^T = 0.1, and changed the correlation length 
£. The correlation time r is much shorter than the quantity ^ACT Xq/ '(fc^T 1 ) 2 so that the correlation between external 
noise / and implicit noise contained in population p is negligible. Here Lq is the typical length scale accompanied 
by the probability distribution. For each parameters, more than 200 trajectories are simulated so that the mean- 
first-passage time well-converges. The time step of the stochastic simulations is around 0.09 which is longer than 
t, but sufficiently short so that the discretization of the equation makes sense. The k dependence of the numerical 
escape rate T is compared with both the theoretical prediction (|17[) and the traditional Kramers formula in Fig. 2. 
We find a very nice agreement between the theoretical and numerical results. In Fig. 3, we also examined the steady 
state distribution for parameters near the maximum of the escape rate of Fig. 2 (r = 0.01, £ = 0.04750, 77 = 1, and 
ksT = 0.1). In this way, the quasi-stable state becomes unstable due to the spatial randomness and leads to the 
thermal renormalisation of both the barrier-height and collisional prefactor in the Kramers escape rate. Noting that 
the temperature is low enough, the normalization factor (1 — At) is positive and the escape rate shows a maximum at 
low but finite temperature. 

In summary, stimulated by the work on the generalized Ornstein Uhlenbeck process in momentum space 0, we 
investigated the diffusion process with the spatially and temporally correlated noise for strong-damping regime. A 
Fokker-Planck equation is derived in a kind of Markovian limit that is characterized by a dimensionless parameter 
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FIG. 2: The k dependence of the escape rate F (|17|) (the solid- line). The broken-line shows the ordinary Kramers escape 
rate(ft = 0). The escape rate is maximized at a finite k. The Landau potential x 4 /4 — x 2 /2 and parameters r = 0.01, r\ = 1, 
and ksT = 0.1 are assumed. The correlation length £ is varied which amounts to the change of ft. 

P st (x) 




FIG. 3: The steady state distribution P st (x) for the Landau potential x 4 /4 — x 2 /2 and parameters are r = 0.01, £ = 0.04750, 
rj = 1, and IcbT = 0.1. The histogram is the result of the stochastic simulation which is composed of 30 samples of 2000 
steps (transient duration is omitted). The solid-line indicates the Eq. (|15[l . and the broken-line is the canonical distribution 
without potential renormalization. 

k = 2rksT/ \/Tr£, 2 rj. The drift term is renormalized by the spatial randomness. Intuitive understanding of this 
phenomenon should be as follows: the spatial randomness yields domains of the spatial coherence of the random noise 
with the typical size £. The systematic force pushes the particle out of the domain of coherence which amounts to 
an extra relaxation of the random noise. Thus the systematic force —U'(x) accompanied by the spatial relaxation 
is equivalently replaced by the weaker force but without the spatial relaxation of the noise. We think that the basic 
scenario is ubiquitous at least qualitatively. The role of the potential renormalization is most clearly seen in the 
escape rate formula. The consequent new escape rate has a maximum at an optimal re, due to a competition between 
the reduction of the prefactor and that of the activation energy. 
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